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Abstract 

We examine the formation of bound states on a generalized nonlinear impurity located at or 
near the beginning (surface) of a linear, tight-binding semi-infinite lattice. Using the formalism of 
lattice Green functions, we obtain in closed form the number of bound states as well as their energies 
and probability profiles, for different nonlinearity parameter values and nonlinearity exponents, at 
different distances from the surface. It is shown that close to the surface, the amount of nonlinearity 
needed to create a bound state or to effect dynamical selftrapping, increases (decreases) depending 
on whether the exponent is smaller (larger) than, approximately, two. 

PACS numbers: 71.55.-i, 73.20.Hb, 03.65.Ge 
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The interplay of nonlinearity and discreteness has received considerable attention 
recently Q], since it plays a vital role in the emergence of a new kind of excitation in extended, 
nonlinear systems with discrete translational invariance, known as intrinsic localized modes 
(ILM). These ILMs are generic to physical systems of interest, such as arrays of nonlinear 
optical waveguides 3, .olecular crystals a, biopol^sQ, a,a yS of JoS epa S oa JU ac tl o ns fl 
and even Bose-Einstein condensates in mag 

Given the strictly local manner in which nonlinearity enters into the effective evolution 
equations in all these cases (see below), one is led to the idea that in the limit of strong 
nonlinearity, one could approximate a typical nonlinear system by a linear one containing a 
small cluster of nonlinear sites, or even a single nonlinear impurity. The system thus sim- 
plified is amenable to exact mathematical treatment, and the influence of other, potentially 
competing effects such as dimensionality, boundary effects, noise, etc., can be more easily 
studied without losing the essential physics. 

For a one-dimensional discrete system in the presence of a single nonlinear impurity, 
located at n = d, the dynamics is given by the well-known discrete nonlinear Schrodinger 
(DNLS) equation: 

.dC„ 



dt 



V(C n+1 + c n _o - x\c n \ p cj n4 , {% = 1) (1) 



where C n is the probability amplitude for finding the excitation on site n, V is the nearest- 
neighbor transfer matrix element, x is the nonlinearity parameter and (3 is the nonlinearity 
exponent. Usually, but not always, /3 — 2 which in a condensed mater context corresponds 
to an underlying harmonic oscillator degree of freedom 'enslaved' to the excitation (electron) 
at the impurity site. When this vibrational impurity is anharmonic in nature, other (3 values 
are possible in principle, with (3 < 2 corresponding to a "hard" vibrational impurity while 
P > 2 corresponds to a "soft" caseQ. 

Bound states for single nonlinear impurities embedded in infinite lattices include 
chaiasQ fl Q, Ca yley tteesQ, tH^Q aad a cuoicQ Q lattices. Now, siace 
the creation of a bound state, or the dynamical selftraping at the impurity site implies the 
localization of energy on a scale of the order of the lattice spacing, one might surmise that, 
by placing the nonlinear impurity at or near the surface of a semi-infinite lattice, the non- 
linearity strength needed to effect localization would decrease, facilitating in this way its 
creation and experimental observation. As an step in that direction, in this work we exam- 
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FIG. 1: A nonlinear impurity near the surface of a one dimensional chain 

ine a simple model consisting of an electron (or excitation) propagating in a semi-infinite, 
linear chain, that contains a single nonlinear impurity at a distance d from the beginning 
('surface') of a semi-infinite chain (Fig. 1), and examine the conditions for the existence of 
bound state (s) and the dynamica 



obtained for the infinite chain 



selftrapping properties, and compare them to the results 



Bound States. We consider Eq.(JIJ for a semi-infinite lattice (n — 0, 1, • • • ) and normalize 
all energies to the half bandwidth of the infinite chain case. The Hamiltonian is given by 

H = 2 £(N< W + M + \n + l)(n\) + l\C d f\d)(d\ (2) 

n=0 

where {\n)} are Wannier states and 7 = x/(2V). The dimensionless Green function G = 
1 f{z - H) can be formally expanded asQ G = G (0) + G^HiG^ + G^ H X G® H X G^ + ■■■, 
where G^ is the unperturbed (7 = 0) Green function and Hi = r )\Cd\ l3 \d){d\. The series 
can be resumed to all orders to yield 

»,\ri 1/3^(0)^(0) 

n _ /^(0) , l\^d\ ^md^dn (o\ 

1-^C^GfJ 

where G mn = (m\G\n). Now, we cannot use Eq.fjSJ) directly since we do not know Cd, but 
we will determine it through an exact selfconsistent procedure: The energy of the bound 
state(s) is obtained form the poles of G mn , i.e., by solving 1 = , y\Cd\^G d ](zi ) )- On the other 
hand, the bound state amplitudes C n are obtained form the residues of G mn at z — z\>. In 
particular, at the impurity site, \Cd\ 2 = Res{Gdd(z)} z =z b = ~G d °] {zb)/G^ (zb) Inserting 
this back into the bound state energy equation leads to 



1_ G d f + \z b ) 
T [-G d { d\z b W 



(4) 



The unperturbed Green function Gmn for the semi-infinite lattice, can be calculated by the 
method of mirror images: Since there is no lattice to the left of n = 0, Gmn should vanish 
identically at n = — 1. Thus, Gmn(z) = G^ n (z) — G^_ n _ 2 (z), where G^ n (z) is the Green 
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function for the infinite lattice, G™ n (z) = sgn(,z)(i/Vz 2 — l)[z — sgn(z)\/z 2 — l]' n m l, where 
sgn(z) = +1(— 1) for z > 0(< 0). Therefore, 

G®(z) = sgn(z)^==iz-sgn(z)vG^l}^ 
y z z — 1 

1 



- sg n(z)^==[z - sgn( 2 ) v / ^l] |n+2+m| . (5) 

From Eq.flHJ) we note the parity property GfJ(-z) = —G^(—z), which implies G^ (— z) = 
G^\—z). This means, according to Eq.fjlJ) that the change 7 — ► —7, reverses the sign 
of Zb. On the other hand, from Eq.(JTJ), it is possible to deduce that the change 7 — > —7 
is equivalent to the change C n — > (— l) n C*. Since we are interested in a localized state, 
where the C n can be chosen as real, we conclude that a change in sign of the nonlinearity 
parameter reverses both the "staggered" character of the bound state and the sign of the 
localized state energy. 

After inserting Eq.(jSJ) into (J3J), the general structure for the number of bound states 
emerges: For any finite distance d from the surface and any positive value of the exponent 
(3, there is a critical amount of nonlinearity 7 below which there is no bound state, and above 
which there are two bound states. For f3 exponents smaller than 2, and as d is increased, one 
of the bound states tends to merge with the band edge, so that in the limit of a very deep 
impurity, there is only a single localized bound state. For /3 > 2, however, as d is increased, 
both bound states remain localized. In the special case of a linear impurity (j3 = 0), there 
is a single bound state provided 7 > l/(d + 1). Thus, in the limit d —> 00 these results 
are consistent with the case single bound state for 

j3 < 2, and for f3 > 2, a critical curve in nonlinearity strength- nonlinearity exponent space, 
separating a region with no bound states from a region with two bound states. At the 
surface (d = 0) the critical curve is given by 

(1 + /3 )(l+/3)/2 

7c = 2^/2 ^ 

In particular, for the DNLS case, 7 C = 3 3 / 2 /4 w 1.3, larger than the value for the infinite 
chain (7 C = 1). Figure 2 shows phase diagrams in 7-/? space showing the number of bound 
states, at different distances d between the impurity and the 'surface' of the system. 

As to the stability of these bound states, it is easy to see from a graphical analysis of 
the structure of Eq.Q that, as nonlinearity 7 is increased one of the bound states becomes 
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FIG. 2: Phase diagram in 7-/? space, showing the number of bound states, for different distances 
impurity-surface (in units of the lattice constant). 
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FIG. 3: Scaled critical nonlinearity for onset of a bound state, as a function of the distance from 
the nonlinear impurity to the 'surface' of the chain. 

more and more localized, while the other becomes more and more delocalized. Since, in the 
limit of high nonlinearity, the effective coupling among sites is negligible, one would expect 
the bound state to become more and more localized. Therefore, the state with smaller 
localization length is stable, the other unstable. This qualitative argument is confirmed by 
the more rigorous procedure of examining the Hamiltonian flow of the system around the 
two fixed points (bound states). 

Figure 3 shows the critical nonlinearity for the onset of a bound state, as a function of 
the distance from the impurity to the lattice 'surface' (n = 0), for different nonlinearity ex- 
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ponents. Significative differences from the infinite lattice case are apparent: As the impurity 
is placed closer and closer to the surface, the critical nonlinearity to create a bound state 
increases or decreases, depending upon whether the nonlinear exponent is below or above, 
approximately two. In particular, for the all-important standard DNLS case, (3 = 2, the 
presence of a surface increases the nonlinearity needed to create a bound state, contrary to 
popular belief that a surface would help localize the electron. 

For a given value of exponent (3 and any inclusion distance d, the bound state probability 
profile \C n \ 2 is given in closed form by \C n \ 2 = A[ Q\ n ^ d \—Q\ n+d+2 \ ] ; where Q = z b — -Jz 2 — 1, 
A = (z b - Q)/{z h + (z b + 2(1 + d)^z\ - l)Q 2(1+d) ) and z b is the solution of Eq.®. Simple 
analysis of this profile shows that \C n \ 2 has always a single hump at n — d. This profile 
is shown in Fig. 4 for the standard DNLS ((3 = 2) and a nonlinearity strength 7 just above 
critical, at four different impurity locations under the surface, its general features are shared 
by other (3 exponents. Below the surface, the probability profile converges quickly to the 
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FIG. 4: Probability profile for the stable bound state at different impurity positions {(3 = 2, 
7 = 1.305). 

infinite lattice d increases past 4. 

Dynamical Seltrapping. We compute numerically the long-time average probability at 
the impurity site, Pd = limT^oo(l/^) Jo dt\Cd\ 2 , for several distances d from the 'surface' 
(n = 0). As initial condition, we use a completely localized excitation on the impurity 
site: C n (0) = 5 n d- Figure 5 shows the critical nonlinearity for selftrapping (Pd > 0) as a 
function of the distance between the nonlinear impurity and the chain surface, for different 
nonlinearity exponent values. In general, the behavior is qualitatively similar to the one 
observed for the onset of a bound state (Fig. 3). In both cases, for a fixed distance, an 
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increase of the nonlinear exponent (3 results on an increase of the nonlinearity threshold for 
selftrapping. The same behavior was observed previously for an impurity in a completely 



infinite chain [ll|. This phenomenon is not hard to explain: Since \C n \ < 1, we see from 
Eq.flTJ) that as (3 > is increased, [C^f will necessarily decrease, which implies that a larger 
7 will be needed to keep the value of the effective impurity strength 7|Cd| /3 . Thus, at a fixed 
impurity-surface distance, a higher (3 implies a higher j c . Another interesting behavior we 
observe from figs. 5 and 3 is that for a fixed nonlinearity exponent, the critical nonlinearity 
depends roughly on whether the exponent is below or above two, approximately: For (3 < 2, 
an increase in the impurity-surface distance d results on a decrease of 7 C , while for (3 > 2, an 
increase in d increases j c . The explanation of this phenomenon seems to rest on the delicate 
balance between kinetic and potential energies. If we assume an electronic bound state *f> 
with localization length A, then on normalization grounds we have |\I/| 2 ~ 1/A. The kinetic 
energy content is AK ~ h 2 /2m\ 2 , while the average potential energy is in magnitude equal 
to AV = j dx V(x)\^(x)\ 2 = J dx 7|\l/(a:)| /3 |\l/(a;)| 2 ~ 7a/A [1+(/3/2)] , where a is of the order 
of the lattice spacing. Thus, 

AV/AK ~ 7 A 1_(/3/2) . (7) 

On the other hand, as the impurity is brought closer to the surface, the wavefunction becomes 
more 'compressed' (Fig.4), i.e., A decreases as d approaches zero. This implies, from Eq.(j2J) 
that for (3 > 2, a decrease in A increases AV^ with respect to AK, which means that less 
nonlinearity is needed to maintain a localized state. On the contrary, if (3 < 2, an decrease 
in A decreases AV with respect to AK and now, more nonlinearity is needed to maintain 
the localized state. 

Completely nonlinear lattice. In the large nonlinearity limit where 7 3> 7 cr , the single 
nonlinear impurity results should approximate those corresponding to a whole nonlinear 
lattice. For the particular case examined in this work, the "extended" problem consists of 
the formation of an intrinsic localized mode (ILM) in a semi-infinite nonlinear latice. Due 
to the presence of a surface, the discrete translational invariance is broken and a natural 
question arises: where will the localized state be formed? Our single nonlinear impurity 
analog can provide an answer. For each impurity position d, the bound state energy can 
be computed as a function of d. The position corresponding to its minimum value will 
correspond to the position of the ILM. Also, the impurity energy and spatial probability 
profile should approximate the ones corresponding to the ILM. Figure 6 shows the impurity 




2 4 6 8 10 
Distance from surface 



FIG. 5: Right: Scaled critical nonlinearity for dynamical selftrapping as a function of the distance 
from the impurity to the chain surface, for several nonlinearity exponents. The empty circles 
shown for (3 = 1 and d = 4 through d = 10 indicate approximate values since the selftrapping is 
not abrupt. 
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FIG. 6: Nonlinear impurity bound state energy as a function of distance impurity-surface, for 
(5 = 2 and 7 = 2 (upper) and 7 = — 2 (lower). 

energies as a function of distance from the lattice surface, for the DNLS case, f3 = 2 and 
7 = ±2. We see that for a positive value of the nonlinearity parameter 7, the preferred 
position is the very surface (d = 0), while for a negative 7, the preferred position is one 
layer below the surface (d = 1). These predictions are indeed confirmed by direct numerical 
computations, where the Hamiltonian corresponding to a semi-infinite nonlinear lattice H = 
(1/2) Y^^=o(\ n ) ( n + 1| + |^ + l)( n l) + 7 Y^=o l^l /3 l n )( n l> i s diagonalized by an iterative 
procedure. For the particular example of fig. 6, the error obtained for the ILM energy is 
about 1%. 



8 



Acknowledgements 

This work was partially supported by Fondecyt grant 1020139. The author is grateful to 
R. Seiringer for useful discussions. 



[1] David K. Cambpell, Sergei Flach, and Yuri S. Kivshar, Physics Today 57, 43 (2004). 

[2] D. N. Christodoulides, R. I. Joseph, Opt. Lett. 13, 794 (1988); Y. S. Kivshar, Opt. Lett. 18, 

1147 (1993); H. S. Eisenberg at al., Phys. Rev. Lett. 81, 3383 (1998); R. Morandotti et al., 

Phys. Rev. Lett. 83, 2726 (1999) and 83, 4756 (1999). 
[3] A. A. Ovchinikov, Sov. Phys. JETP 30, 147 (1970); R. Bruinsma et al., Phys. Rev. Lett. 57, 

1773 (1986); A. S. Davydov, Solitons in Molecular Systems, E. S. Kryachko, trans., Kluwer, 

Hingham, Mass. (1985). 

[4] A. Xie et al., Phys. Rev. Lett. 84, 5435 (2000); T. Dauxois, M. Peyrard, in Nonlinear Exci- 
tations in Biomolecules, M. Peyrard, ed., Springer, New York (1995), p. 127; J. C. Eilbeck 
at al., Physica D 16,318 (1985); A.C. Scott, Phys. Rep. 217, 1 (1992); A. Scott, Nonlinear 
Science: Emergence and Dynamics of Coherent Structures, 2nd ed., Oxford U. Press, New 
York (2003); G. P. Tsironis at al., Europhys. Lett. 57, 697 (2002); S. F. Mingaleev at al., 
Europhys. Lett. 59, 403 (2002). 

[5] L. M. Floria et al., Europhys. Lett. 36, 539 (1996); E. Trias et al., Phys. Rev. Lett. 84, 741 
(2000); P. Binder et al., Phys. Rev. Lett. 84, 745 (2000); A. Ustinov, Chaos 13, 716 (2003). 

[6] A. Trombettoni and A. Smerzi, Phys. Rev. Lett. 86, 2353 (2001); J. Phys. B 34, 4711 (2001); 
A. Smerzi, A. Trombettoni, P.G. Kevrekidis, and A.R. Bishop, Phys. Rev. Lett. 89, 170402 
(2002). 

[7] M. I. Molina, Mod. Phys. Lett. B, 17, 1 (2003) and unpublished. 

[8] O. M. Braun and Y. S. Kivshar, Phys. Rep. 306, 1 (1998); A. A. Sukhorukov, Y. S. Kivshar 
and O. Bang, Phys. Rev. E 60, 41 (1999); R. A. Vicencio, M. I. Molina and Y. S. Kivshar, 
Phys. Rev. E, in press. 
[9] D. Cheng, M. I. Molina and G. P. Tsironis, J. Phys.: Condens. Matter 5, 8689 (1993). 
[10] M. I. Molina and G. P. Tsironis, Phys. Rev. B 47, 15330 (1993); 
[11] G. P. Tsironis, M. I. Molina and D. Hennig, Phys. Rev. E 50, 2365 (1994). 



9 



[12] M.I. Molina in Topics in Theoretical Physics, edited by V.C. Aguilera-Navarro, D. Galletti, 
B.M. Pimentel and L. Tomio, (IFT, Sao Paulo, 1996); M.I. Molina, Phys. Rev. B 60, 2276 
(1999). 

[13] M.I. Molina, Phys. Rev. B 60, 2276 (1999); K. M. Ng, Y.Y. Yiu and P.M. Hui, Solid State 

Commun. 95, 801 (1995). 
[14] B. C. Gupta and S. B. Lee , Phys. Rev. B 63, 144302 (2001). 
[15] Y. Y, Yiu, K. M. Ng and P. M. Hui, Phys. Lett. A 200, 325 (1995) 
[16] C. A. Bustamante and M. I. Molina, Phys. Rev. B 62, 15287 (2000). 
[17] M. I. Molina, work in progress. 

[18] E.N. Economou, Green's Functions in Quantum Physics, Springer Series in Solid State Physics, 

Vol.7 (Springer- Verlag, Berlin, 1979). 
[19] M. I. Molina, Mod. Phys. Lett. B 17, 1 (2003). 



10 



